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SUMMARY 

The Dynamic Engine ANalysIs program, DEAN, Is a FORTRAN code Implemented 
on the IBM/370 mainframe at NASA Lewis Research Center for digital simulation 
of turbofan engine dynamics. DEAN Is an Interactive program which allows the 
user to simulate engine subsystems as well as a full engine systems with rela- 
tive ease. The non-linear first order ordinary differential equations which 
define the engine model may be solved by one of four Integration schemes, a 
second order Runge-Rutta, a fourth order Runge-Kutta, an Adams Predictor- 
Corrector, or Gear's method for stiff systems. The numerical data generated 
by the model equations are displayed at specified Intervals between which the 
user may choose to modify various parameters affecting the model equations and 
transient execution. Following the transient run, versatile graphics capa- 
bilities allow close examination of the data. 

DEAN'S modeling procedure and capabilities are demonstrated In this paper 
by generating a model of a simple compressor rig. 


INTRODUCTION 

System simulation Is a powerful tool In the aerospace Industry. These 
simulations can be used In conjunction with testing of a specific system In an 
attempt to anticipate possible problems ana reduce risk. Also simulations can 
be utilized In subsystem development, specifically In the design and develop- 
ment of controls for turbofan engines. 

These simulations must have several attributes In order to be useful 
tools. They must be reasonably accurate and computationally efficient. 

Currently there are several programs In use that simulate turbofan 
engines. NNEP (ref. 1) Is a generalized code for steady-state analysis of 
arbitrary engine configurations. Another code, DYNGEN (ref. 2) simulates the 
dynamics of specific engine configurations. A complex hybrid (analog/digital) 
computer simulation, HYDES (ref. 3) simulates dual shaft turbofan engines with 
up to three air streams. These simulations are limited as to their complexity 
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because the hybrid computer has hardware and memory limitations. D1GTEM 
(ref. 4), a recently written code allows steady-state as well as transient 
calculations to be performed. DIGTEM has the flexibility to model various 
engine configurations which are subsets of a two-spool, two stream turbofan. 
DIGTtM's Interactive capability Is limited and turbofan engine components can- 
not be developed separately. None of the above packages have the extensive 
flexible graphics provided In DEAN. 

The Dynamic Engine ANalysIs program, DEAN has been written specifically 
to model two-spool, two-stream turbofan engines and turbofan engine components. 
However, DEAN can be easily reconfigured to model other dynamic systems such 
as turboshaft engines. Internal combustion engines, or electrical power 
systems. User Interaction, modular engine component routines, various Inte- 
gration methods, and extensive graphics have been Incorporated Into the soft- 
ware to enhance the simulation capability of DEAN. 

Implemented under the TSS operating system on an IBM/370 mainframe com- 
puter, DEAN was written with the Intent of simulating dynamic systems Inter- 
actively. The program utilizes the IBM/370 ' s Interactive debugging 
capabilities, graphics libraries, and the IMSL FORTRAN library which allows 
the user to develop models on-line In the sense of finding modeling errors and 
Initiating transients. The Interactive capability also allows the user to 
analyze system response to external Inputs such as perturbations In fuel flow, 
nozzle area, and bleed flows. Temperatures, pressures, flows, and other model- 
ing parameters otherwise known as state variables are displayed continuously 
during the simulation. Also, part of DEAN'S user Interaction Includes Its user 
friendliness which Is achieved by using menus with logical default selections 
during program execution. 

There Is a natural division of components In a turbofan engine, l.e., 
Inlet, fan, compressor, combustor, high pressure turbine, low pressure turbine, 
afterburner, nozzle, bypass duct, and engine control. DEAN Is designed to 
simulate a turbofan engine with a set of modular routines, one for each engine 
component. The user may use the provided component models for his simulation 
or he may choose to provlde/develcp his own model for any or all of the 
components . 

Generally the models are based on first principles, l.e., conservation of 
mass, momentum, and energy, and can be formulated as a system of first order 
nonlinear differential equations. DEAN uses these equations and one of several 
numerical Integration techniques to approximate the system response. The 
Integration technique may be selected by the user which best approximates the 
system being simulated. The appropriate use of the available Integration 
schemes will be discussed later. 

DEAN supports an extensive graphics capability for displaying results of 
a simulation. This capability Includes the ability to plot any parameter, 
Including time, versus any other parameter, superposition of performance maps 
on a plot, calculation of new parameters based on existing parameters, the 
option of having one or two plots per page, extensive labeling Including date, 
time, and confidentiality of data, ability to plot the data for time sub- 
intervals, and the capability to save the data for later analysis. The plots 
generated by DEAN are of report quality. 
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There are several support packages required to model turbofan engines. 
Generally the engine performance maps are bivariate In nature, l.e., compressor 
pressure ratio Is a function of corrected speed and corrected flow. Also var- 
ious thermodynamic properties are required such as C p (air and fuel), C v , 
and y, the ratio of specific heats, all of which are single valued functions. 
Therefore subroutines have been supplied to Interpolate single and bivariate 
maps. Also there are subroutines for reading free formated datasets and a sub- 
routine Is Included for simple proportlonal/lntegral control. 


DEAN STRUCTURE 

DEAN was designed to be a modular program. Subroutines were written to 
do specific tasks such as prompting the user for a steady-state point or run- 
ning a transient simulation. This manner yields several advantages. DEAN can 
be expanded and modified easily because only specific modules may need minor 
changes. Also, the modularity In the program expedited the software develop- 
ment phase. Further, the programming flow of DEAN Is more easily understood. 
Thus DEAN Is a highly structured program which allows the user to easily Iden- 
tify and modify components. 

DEAN consists of 46 subroutines that are called by one main program. 

Since there are so many subroutines, only the more Important ones will be dis- 
cussed here. Figure 1 shows a simplified block diagram of the subroutine 
Interaction In DEAN. 



Figure 1. - Simplified diagram of Duan subroutines. 
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The main program (DEAN) Is responsible for controlling the entire program. 
It presents the user with the following options: (1) system definition, (2) 

map definition, (3) steady-state operating point calculations, (4) transient 
simulation, (5) data analysis (graphics), (6) termination of DEAN - and prompts 
for a response. When the user has entered an appropriate response, the program 
calls the subroutine associated with that response. 

Subroutine DEFSYS Is prompts the user for the components of the system to 
be simulated. By using a menu of components (table 1.), DEFSYS allows the user 
to define a system for simulation. The order of these components corresponds 
to the configuration of a turbofan engine. If the user desires to simulate a 
system other than a turbofan engine, he must redesignate the components to 
correspond to his system. For example, when simulating a turboshaft engine, 
the low pressure turbine component would contain a model for torque loading. 
With minimal software modifications to various menus In DEAN, the system con- 
figuration may be formally adapted. The user will also be prompted for the 
required model conditions such as Inlet temperature, pressure, and fuel flow 
If the combustor and/or the afterburner are being simulated without a fuel 
control . 

From the components chosen by the user, subroutine MAPSET decides which 
bivariate maps may be needed for the simulation and prompts the user to execute 
one of two tasks related to Initialization of these maps. The rough map data, 
supplied by the user, resides In specified format In a dataset. When beginning 
a simulation MAPSET processes the rough data Into a uniform grid for the bivar- 
iate Interpolation routines and, upon request, rewrites the data over the 
original dataset. The processing of rough data needs to be done only once 
because MAPSEI can determine when the data In the dataset Is rough data or the 
pre-processed grid. DEAN Is capable of processing eight bivariate maps which 
have been designated as shown In table 2. The user may wish to change the 
designated maps and, may do so as long as no more than eight bivariate maps 
are specified. Also, upon request, MAPSEI will graphically display the 
processed map data as each map Is processed. This allows the user to verify 
the maps specified by the data he has supplied. 


TABLE 1. - LIST OF COMPONENTS 
SIMULATED BY DEAN 


Component number 

Turbofan engine component 

1 

Inlet 

2 

Fan 

3 

Compressor 

4 

Combustor 

5 

High pressure turbine 

6 

Low pressure turbine 

7 

Afterburner 

8 

Nozzle 

9 

Bypass 

10 

Rotor dynamics 

11 

Main fuel control 
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TABLE 2. - AVAILABLE BIVARIATE MAPS 
SUPPORTED BY DEAN 


Map number 

Description 

1 

Fan pressure ratio 

2 

Fan temperature ratio 

3 

Compressor pressure ratio 

4 

Compressor temperature ratio 

5 

High pressure turbine pressure ratio 

6 

High pressure turbine enthalpy 

7 

Low pressure turbine pressure ratio 

8 

Low pressure turbine enthalpy 


Subroutine SSOPPT Is responsible for providing the simulation with an 
Initial starting point or steady-state point. There are two mechanisms the 
user may choose to determine such a point. Generally this point consists of a 
vector of state variables that represent the pressure, temperatures, flows, 
shaft speeds, etc. of the system. There Is one state variable for each differ- 
ential equation In the system model as well as any parameters calculated by 
empirical or theoretical equations such as choked flow for example. The first 
mechanism Is to Input the states directly from the keyboard. DEAN will prompt 
the user for the states, component by component as specified by the system 
definition portion of DEAN. An option Is provided to use this Input as an Ini- 
tial guess In conjunction with the governing equations and Iterate to a steady- 
state point by a modified Newton- Raphsrn technique. An alternate method for 
determining a steady point Is to use the InUlal guess as a starting point for 
a transient and run the transient until the dynamics settle out. Once a steady 
state point has been calculated, It may be stored with the system configuration 
for future reference. This allows the user to select a saved steady-state 
point rather than entering the states each time a simulation Is run. 

The subroutine TRNSNT Is responsible for running the transient simulation. 
TRNSN1 prompts the user for one of four Integration schemes and, after an 
acceptable Input, requests more specific Information such as Initial time step, 
length of transient time, frequency at which states should be displayed at the 
terminal, and also whether the user wishes to Interact directly with the sim- 
ulation. Data generated by the transient simulation Is stored at specified 
Intervals for graphical analysis at a later time. 

The subroutine PLOT Is responsible for generating graphics output from 
DEAN. This subroutine supports various options highlighted by table 3. These 
options are discussed briefly later In the paper. 

The last group of Important subroutines are listed In table 4. These sub- 
routines contain all of the modeling equations required by the simulation. 
Included In these subroutines are data for the component geometry, data for 
single variable functions, code to store data for graphical analysis, bivariate 
map calls, evaluation of specific thermodynamic properties such as 
Cp = f (temperature) , etc. These subroutines can be written by the user who 
wishes to develop an analytical model, as opposed to using those provided. 
Currently DEAN Is set up to Integrate up to one hundred state variables. 
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TABLE 3. - GRAPHICS PACKAGE OPTIONS 
FROM SUBROUTINE PLOT 


Option 

Description 

1 

Read plot labels from labels dataset 

2 

Read plot data from dataset 

3 

Save plot data on dataset 

4 

Define or modify plot parameters 

5 

Display a menu of available plot parameters 

6 

Change the time Interval of the data plotted 

7 

Automatically make a hardcopy of each plot 

8 

Generate a standard cartesian plot 

9 

Superimpose a map on a standard cartesian plot 

10 

Put data classification labels on all plots 


TABLE 4. - LIST OF DEAN'S 
COMPONENT SUBROUTINES 


Subroutine 

Description 

EQUNS 

Subroutine for "global" equations 

INLET 

Inlet duct model 

FAN 

Fan model 

COMP 

Compressor model 

BURNER 

Combustor model 

HPTURB 

High pressure turbine model 

LPTURB 

Low pressure turbine model 

AFTBRN 

Afterburner model 

NOZZLE 

Nozzle model 

BYPASS 

Bypass duct model 

MFC 

Engine control model 

ROTDYN 

Dynamic rotor balance equations 


It may be noted, that while the component subroutine names correspond to 
those of a turbofan engine, they may also be used to model components of vir- 
tually any dynamic system. 


INTERACTIVE ASPECTS OF DEAN 

The simulation program, DEAN, has five major areas of user Interaction. 
The first area Involves defining the system components (listed In table T). 
DEAN will prompt the user for the types of components required for the Simula 
tlon, and thus subsystems of the turbofan engine may be modeled, l.e., compres 
sor and nozzle. Based on the component models used, DEAN will prompt the user 
for specific Information such as Inlet pressure and temperature of the model 
or fuel flow for a combustor. 

The second area of Interaction Is map grid definition for bivariate maps 
(z = f[x,y]). DEAN will prompt the user for Information concerning the deflnl 
tlon of the grid. The program can presently handle as many as eight bivariate 
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maps as listed In table 2. Such maps are usually required for fan, compressor, 
high pressure turbine, and/or low pressure turbine performance data. 

The next area of Interaction concerns the Initial values of the states for 
the simulation. The user must Initially enter values for the states which give 
the Integration schemes a starting point. These values, which usually corre- 
spond to a steady point are entered from the keyboard. The user may also 
choose to select a previously entered steady point. 

The most useful area of Interaction Is during the actual system simula- 
tion. As the transient proceeds, the program will save data on specified 
parameters (up to 100) with a maximum of 1000 points per parameter. The method 
of saving the data Is an Integral part of the component subroutines. During 
the transient the user may change various parameters based on the simulation's 
output. 

The final area of Interaction Is In data analysis. DEAN supports a menu 
driven graphics package (table 4) with capabilities for basic mathematical 
manipulation of the transient data. 


MODULAR ENGINE COMPONENTS 

A turbofan engine Is generally comprised of those components listed In 
table 1. In the Dynamic Engine ANalysIs program there Is a subroutine allo- 
cated for each of these components. The subroutines consist of FORTRAN code 
for evaluating derivatives and storing data for analysis. Also, there Is code 
to generate parameters, perform map lookups for single and bivariate maps, and 
evaluate thermodynamic properties for specific components. These component 
subroutines are always called by DEAN In the order shown In table 1. 

This modular simulation structure allows the user to develop turbofan 
engine systems and subsystems, component by component. This method Is useful 
because It allows easier debugging of both coding and theoretical problems. 
Currently, the only restriction Is that any subsystem must have Its components 
ordered as per table 1. For example, a turbofan core consisting of a compres- 
sor, a combustor, and a turbine In line may be simulated. However, a change 
In component arrangement or system topology would require minor code changes 
and possible renaming or addition of component subroutines. 


INTEGRATION SCHEMES 

The mathematical model of the system Is represented generally by ordinary 
first order nonlinear differential equations. There are several methods which 
can be used to Integrate these equations. DEAN supports the following Inte- 
gration schemes: (1) second order Runge-Kutta, (2) fourth order Runge-Kutta, 

(3) Adams predictor-corrector, and (4) Gear's method for stiff systems. Each 
of these Integration schemes has Its own appropriate applications. 

The second and fourth order Runge-Kutta schemes (ref. 5) execute simi- 
larly. The second order scheme requires fewer function evaluations, and thus, 
less CPU time. However, the fourth order scheme Is more accurate. In general, 
these explicit methods should be used when model frequencies and accuracy 
requirements are high. 
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The Adams predictor-corrector Integration (ref. 6) has several appealing 
attributes. It Is a variable order (up to 12), variable time step method. The 
algorithm Is a multlstep process In which there Is a predictor for extrapolat- 
ing to the solution and a corrector for Interpolating a solution based on the 
extrapolation. The Interpolation procedure Involves Iteration and, therefore, 
requires more computation time than the one-step Runge-Kutta (fourth order 
method) If there are more than two Iterations In the corrector process. 

The Gear's algorithm for Integrating differential equations (ref. 7) was 
developed for systems that have a wide range of frequencies (stiff systems). 

A turbofan engine can be considered as such a system. The Gear's method Is a 
backward difference scheme capable of variable order (up to five) and variable 
time step. It Is useful for simulations that require a substantial amount of 
computer time for function evaluations. 


GRAPHICS CAPABILITIES 

The last of DEAN's main menu options Is the Data Analysis package. This 
set of subroutines allows the user to manipulate, as well as plot, data gen- 
erated during a transient. The main subroutine which controls execution of 
this set of code Is subroutine PLOT. PLOT contains a menu of ten options for 
the user to choose from. 

The ten PLOT options allow the user to: (1) read user specified plot 

labels corresponding to the parameters saved during a transient, (2) read plot 
data stored In a dataset from a previous transient, (3) write transient data 
Into a dataset for future reference, (4) perform basic mathematical manipula- 
tion of the transient data and/or modify plot title and labels, (5) display a 
list of the plot labels, (6) change time Interval over which the data Is plot- 
ted, (7) automatically hardcopy plots (assuming hardcopy capability exists), 

(8) produce cartesian plots of specified parameters, (9) plot specified para- 
meters with a bivariate map (corresponding to table 2) super-imposed on the 
data, (10) write security classification labels on plots of the transient data. 


CONCLUDING REMARKS 

DEAN Is a versatile simulation program written In FORTRAN IV. Designed 
for Interactive dynamic simulation of turbofan engine systems and subsystem , 
DEAN has successfully simulated a number of component models, such as the com- 
pressor rig example In the Appendix and a full turbofan engine. Further, the 
Interactive graphics capabilities have proven to be useful for simulation 
development and analysis. 

DEAN Is supported by Interactive debugging, the IMSL library (9th 
edition), and a FORTRAN accessible graphics package. 
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APPENDIX A 


COMPRESSOR RIG MODEL EXAMPLE 


The Dynamic Engine ANalysIs program can best be understood by actually 
developing a simulation of a system. Here, a simple simulation of a compressor 
rig Is used to demonstrate the capabilities of DEAN. Also, the model develop- 
ment will not overshadow the use of DEAN as a simulation program. 

First the model Is defined. Let us assume that the compressor rig can be 
modeled as a constant speed, one-dimensional, multi-lumped volume, Incompress- 
ible system with three major components: (1) Inlet volume, (2) compressor, and 

(3) turbine nozzle. Note that these three elements correspond to those In 
table 1. The Inlet will be modeled as a single lumped volume and the comores- 
sor will also be represented as a single volume employing the actuator disk 
concept suggested by Greltzer (ref. 8) for compressor stall. The turbine 
nozzle will be modeled by using a choked flow relation between temperature and 
pressure. 

Generally three parameters provide communication between components: here 

(l)pressure (P), (2) mass flow (ft), and (3) temperature (T) are used. 

These parameters are calculated by Integrating the non-linear first order dif- 
ferential equations that result from the conservation of mass, momentum, and 
energy and the equation of state In differential form. Figure 2 shows a typi- 
cal lumped volume with the various parameter designations. The general form 
of these equations for such a lumped volume analysis Is as follows, where p 
Is density, t Is time, V 0 i Is volume, R Is gas constant of air, L Is the 
volume length, A Is cross sectional area, F ex ^ denotes sum of frictional 
and pumping forces, Q denotes the net heat transfer to the volume, W denotes 
net work done by the volume, C v Is the constant volume specific heat of the 
fluid, and y the ratio of specific heats. 
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figure 2. * Typical lumped volume diagram. 
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figure 3. - System block diagram. 


TABLE 5. - STEPS REQUIRED TO DEVELOP A SIMULATION 


Step number 

Instructions 

T 

Define the system to be modeled 

2 

Divide the system Into components 

3 

Write dynamic equations for each component 

4 

Write FORTRAN component subroutines 

5 

Create datasets to be used by simulation 

6 

Access and execute DEAN 
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From these equations, models of an Inlet and compressor can be developed. 

Table 5 shows the general steps Involved In developing the simulation. A block 
diagram of the system Is shown In figure 3. Using the block diagram and the 
general conservation equations the Inlet equations may be written as follows: 


pressure (state) 



V ol.l 






mass flow (momentum) 


dW 

dt * L, 


P 1 A 1 - 


P 2 A 2 " 


r w: 


terr^rature (energy) 


dT, 


RT, 


4t Vo). 2 


S, ( -r T , - Tj) - i c (-,l 2 - t 2 ) 


where F ex t Is modeled as F ex t » -RW 2 (R = constant to give appro- 
priate pressure drop). Also, for the Inlet, P and T are constant Inlet 
conditions. The compressor Is modeled In much the same way. The equations for 
the compressor are as follows: 

d p . RT P dT 

pressure (state) ^ - rTZ <* c - V * t it 


0* c 9 C | 

mass flow (momentum) * — |P ? A 2 - P 3 A 3 ♦ F pxt 


F ext * P 3,map A 3 P 2,map A 2 


temperature (energy) 



RT, 


P V 

r 3 ol ,2 



V - *„<yT 3 - Ij) 


Y ^ T 3,map ” ^.map^j 

where the values of Pj>,map and p 3,map are the pressures corresponding 
to a steady state compressor performance map In terms of pressure ratio (pres- 
sure ratio Is a function of corrected flow and corrected speed) and T 2 
and T 3,map are temperatures corresponding to another steady-state perform- 
ance map In terms of temperature ratio (temperature ratio Is a function of 
corrected flow and corrected speed.) The turbine nozzle will be modeled by 
using the following choked flow nozzle equation: 
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The governing equations are now coded Into separate subroutines for each 
component. Once these subroutines are written for all of the component mcdels 
the user must collect all the data required for the model. For the compressor 
example there are two data requirements; (1) compressor map data and (2) geom 
etry data. 

Once all of the preliminary work Is done for the simulation DEAN may be 
run. A typical run Is similar to that discussed In the INTERACTIVE ASPECTS OF 
DEAN section. 

ror this particular example, the nozzle exit area Is ramped closed during 
the transient while compressor shaft speed Is held constant as shown In 
figure 4. The Adam's method was chosen as the Integration scheme with an 
Initial time step of 0.0001. The length of the transient was 0.50 sec which 
required 6.63 sec of CPU to Integrate six states. This disturbance causes a 
compressor Instability called surge. Figure 5 shows graphically the Instabil- 
ity superimposed on the appropriate performance map. Figures 6 through 8 show 
time traces of Inlet and compressor exit flow, temperature, and pressure, 
respectively. 
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Figure 5. - Plot of compressor instability. 


NON-DIMENTIONAL NON-DIMENTIONAL NON-DIMENSIONAL NON-DIMENSIONAL 

COMPRESSOR EXIT NOZZLE WEIGHT FLOW COMPRESSOR WEIGHT INLET WEIGHT FLOW 

TEMPERATURE FLOW 


ORIGINAL FAClw Sis' 

OF POOR QUALITY 


<•» D E A N -- DYNAMIC ENGINE ANALYSIS PROGRAM *• 



Hi 


TIME 

DEAN: APR. 26. 1985 NASA LEWIS RESEARCH CENTER 
CLEVELAND. OHIO 

Figure 6. - Time traces of weight flow. 
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Figure 7. - Time traces of weight flow and temperature. 
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Figure 8. - Time traces of pressure. 
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